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Abstract 

The conformal Einstein equations for a tracefree (radiation) perfect fluid are derived in 
terms of the Levi-Civita connection of a conformally rescaled metric. These equations are 
used to provide a non-linear stability result for de Sitter-like tracefree (radiation) perfect fluid 
Friedman-Lemaitre-Robertson-Walker cosmological models. The solutions thus obtained ex- 
ist globally towards the future and are future geodesically complete. 
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ON . 

\0 ■ 1 Introduction 

The conformal Einstein field equations have proven a powerful tool to analyse the stability and 
the global properties of vacuum, electro- vacuum and Yang-Mills-electro- vacuum spacetimes — see 
e.g. [S] Uni [HI [m Uni I2D] • By contrast, to the best of our knowledge, there has been no attempt 
to make use of conformal methods to analyse similar issues in spacetimes whose matter content 
is given by a perfect fluid. In this article we make a first step in this direction. We discuss 
the stability and the global properties of a class of cosmological spacetimes having as a source a 
perfect fluid with tracefree energy-momentum tensor. The solutions we construct are non-linear 
' perturbations of a Friedman-Lemaitre-Robertson- Walker (FLRW) reference spacetime. 

The present analysis is to be regarded as a first step in the development of conformal methods 
for the discussion of cosmological models whose matter content is described by a perfect fluid. 
Hence, we restrict our attention to the simplest case from the point of view of conformal methods: 
perturbations of a traceless prefect fluid cosmological model with compact spatial sections of 
positive constant curvature. Generalisation of our analysis to more general background solutions 
and equations of state will be discussed elsewhere. 

The problem of the non-linear stability of FLRW cosmologies and the exponential decay of 
perturbations is considered in |23j . In that reference, a frame formulation of the Einstein-perfect 
fluid system [12J is used to obtain a suitable symmetric hyperbolic evolution system for which the 
Kreiss-Lorenz theory can be readily applied — see [T7]. The results obtained hold for a large class 
of equations of state, but not very stiff ones — like the pure radiation case discussed in the present 
article. More recently, the problem of the non-linear stability of the irrotational Euler-Einstein 
system for de Sitter-like spacetimes has been analysed in [24J. This analysis shows that FLRW 
background solutions with pressure p and density p related by a barotropic equations of state 
of the form p = jp with 1 < 7 < | are future asymptotically stable under small irrotational 
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perturbations. An extension of this analysis to the case of fluids with non-zero vorticity has been 
given in [26J. 

It is notable that the case of a pure radiation perfect fluid cannot be covered by the analysis 
of p51 [Ml By contrast, from the point of view of conformal methods, the pure radiation 
perfect fluid case turns out to be one of the simplest scenarios to be considered. Finally, it 
should be mentioned that conformal methods have been used to pose an initial value problem 
for the Einstein-Euler system at the Big Bang for a class of cosmological models with isotropic 
singularities — see [T]. The methods used in that work do not allow, however, to obtain global 
existence assertions towards the future. 

Our main result can be stated as follows: 

Theorem. Suppose one is given Cauchy initial data for the Einstein-Euler system with a de 
Sitter-like cosmological constant and equation of state for pure radiation. If the initial data is 
sufficiently close to data for a FLRW cosmological model with the same equation of state, value 
of the cosmological constant and spatial curvature k = \, then the development exists globally 
towards the future, is future geodesically complete and remains close to the FLRW solution. 

A detailed and technically precise version of this result is given in Theorem [51 

Remark 1. Similar future global existence and stability results can be obtained using the meth- 
ods of this article for a FLRW background solution with pure radiation equation of state, de 
Sitter- like or vanishing cosmological constant. A, and k = 0, —1. These models expand indef- 
initely towards the future, and remarkably, their scale factor can be computed explicitly — see 
|14| . In the cases with A = 0, minor technical modifications need to be introduced to account for 
a null conformal boundary. The stability of these models will be discussed elsewhere by means of 
different (conformal) methods. 

Remark 2. The restriction of our analysis to the case of perfect fluids with traceless energy- 
momentum tensor is technical: in this case the equation of conservation of energy momentum 
transforms homogeneously under conformal transformations. In the case of perfect fluids with an 
energy-momentum tensor with non- vanishing trace a regularisation of the rescaled equations of 
motions must be carried out. The analysis for the wave equation in [21 115| may be a guide for 
this type of generalisation of our analysis. 

Structure of the article 

The article is organised as follows: Section[2]provides a summary of the tensorial conventions to be 
used in the present article. Furthermore, in Subsection l2.2l a discussion of the procedure of how to 
coordinatise and introduce frame fields of the 3-sphere, S"^ is provided. Section [3l provides general 
remarks concerning perfect fiuid cosmological models and a summary of the properties of the 
background solutions required in our subsequent analysis. These are summarised in Proposition 
[Tl Section [4] gives a brief summary of the conformal Einstein field equations with matter. Section 
[51 provides a discussion of the Euler equations in the context of the conformal field equations. In 
Section [6l we discuss gauge considerations and the procedure leading to a hyperbolic reduction of 
the conformal field equations. The keys steps in this procedure have been discussed extensively 
elsewhere, so that this discussion is kept to a minimum. In particular. Subsection 16.21 provides a 
summary of the structural properties of the conformal evolution equations while Subsection 16.31 
analyses the issue of the propagation of the constraints. Section [71 casts the FLRW background 
as a solution of the conformal field equations of Section [4l and analyses some of its properties. 
Finally, Section[8lis concerned with our main result — the existence and stability result for perfect 
fluid cosmologies with a de Sitter-like cosmological constant as given in Theorem [51 
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2 Notation and conventions 



2.1 Index and curvature conventions 

Throughout this article we work with a spacetime (A^,^^^), where g^u, (m, = 0,1,2,3) is a 
Lorentzian metric with signature (+,—,—,—). We will denote by V the Levi-Civita connection 
of g^Jl,u — that is, the unique torsion- free connection that preserves the metric g^u- In the sequel, 
Rp,i^\p, Rfiu and R will denote, respectively, the Riemann curvature tensor, the Ricci tensor and 
the Ricci scalar of the Levi-Civita connection V. The conventions for the curvature used in this 
article are such that 

i?^.Apr= (VaVp-VpVa)^^, i?M- = i?\AM, R^R^.r"- (1) 

As a consequence of our signature conventions, then A < corresponds to de Sitter-like values of 
the cosmological constant, while A > corresponds to anti-de Sitter-like values. While /i, i', . . . 
denote spacetime tensorial indices, a, (3, ■ ■ ■ denote spatial tensorial ones. Most of our discussion 
will be based on a frame formalism in which i, j, . . . denote spacetime indices ranging 0, . . . , 3. 
Similarly, a, 6, ... will denote spatial indices ranging 1, 2, 3. Spinorial expressions and arguments 
will be used routinely, and we will follow the conventions of [21J. Consequently, the indices 
A, B, . . . will be spinorial ones. 



2.2 Coordinates and vector fields on the 3-sphere 

The present analysis will be concerned with spacetimes which are conformal to manifolds with 
topology / X S'^ where / is an open interval on M. In what follows, the manifold S"^ will always 
be thought of as the following submanifold of R"* : 



The restrictions of the functions x-^, A — 1, 2, 3, 4 on M"* to S'' will again be denoted by x^. The 
vector fields 

Ci=x^d4-x'^di+x^d3-x^d2, (2a) 
C2 = x^d3-x^di + x'^d2-x^d4, (2b) 
C3 = x^d2-x^di + x^d4-x'^d3, (2c) 

on are tangent to S'^. In the sequel, they will always be considered as vectors on S'^. The 
vector fields {c^} = {ci, C2, C3} constitute a globally defined frame on S'' which is orthonormal 
with respect to the standard metric of Moreover, the frame {ci, C2, C3} can be completed 
with a vector cq which is orthonormal to the standard metric on / X {Cs} = {cq, Ci, C2, C3}. 

Let {M.,gfiw) be a spacetime such that the manifold M is diffeomorphic to R x S"^. A map <i> 
defined on an open subset U G A4 will be said to be a cylinder map if it maps U diffeomorphically 
onto a set / X §^ , such that the sets <I>~"'^({r}xS^) are spacelike Cauchy hypersurfaces of and the 
curves / 9 r — )■ $~^(t,p) C A^, p G §^ are timelike with respect to the metric g^i,. The cylinder 
map will be used to pull-back to ti the coordinates (r, a;-^) = {t,x) in / x S'^. Furthermore, one 
can use $ to pull-back to U the frame fields Cg_ defined in the previous paragraph. For simplicity 
of notation, such pull-back will be denoted again by c^. 



3 General remarks about FLRW cosmological models 

A cosmological model [Ai , g , u^) is a representation of the universe at a particular averaging 
scale. It is defined by a Lorentzian metric g^^ on the manifold M and by a family of fundamental 
observers whose congruence of worldlines is represented by the timelike 4-velocity u'' — usually 
taken to be the matter 4-velocity. It is usually assumed that this congruence is expanding at 
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some time. These assumptions together with a specification of the matter content are used to 
determine the dynamics of the universe. In what follows, it will be assumed that the interaction 
between geometry and matter is described by the Einstein field equations 

Rfj-u ~ \RgtJ.iy + •^.9^11' — (3) 
and the energy-momentum conservation equation 

- 0. (4) 

As already mentioned, the conventions for the cosmological constant A used in the present article 
are such that in vacuum, the case A < describes a de Sitter-like spacetime, while the case A > 
corresponds to an anti-de Sitter-like one. 

Our discussion will be concerned with energy-momentum tensors of perfect fiuids for which 

r^i/ = {p + p)Uf^u^ ~ pgi_,y, 

where p, p and -0'' denote, respectively, the density, pressure and 4-velocity of the cosmological 
fluid. The fluid 4-velocity u'^ is timelike and satisfles the normalisation condition u^u^ = 1. 

The background solution whose non-linear stability will be considered in the present article 
belongs to the family of so-called Friedman-Lemaitre-Robertson- Walker (FLRW) cosmological 
models. The FLRW models are homogeneous and isotropic. Their line element is usually given 
in the form 

= -g^Ax^Ax" = dt^ "'^^^^ , (dr^ + r'^AO'^ + gj^^2 g,^^2^ ^ 

(l + ifcr2)- 

where a[t) is the so-called scale factor. This metric automatically defines a perfect fluid energy- 
momentum tensor. When fc = the spatial sections are flat, if A: < the spatial sections have 
negative curvature, while if fc > 0, the spatial sections have positive curvature. The present 
analysis is concerned with FLRW cosmologies with spatial sections of positive curvature (k = 1) 
for which coordinates can be introduced such that: 

~g^ = At'-a''{t)da^, (6) 

with 

dcr^ = di/"^ -I- sin^ ipdO'^ + sin^ ip sin^ Odip^ , 

the standard line element of S'^ in polar coordinates. If the cosmological fluid satisfles the 
barotropic equation of state p = (7 — l)/5, where 1 < 7 < 2 is a constant, then the evolution of 
the scale factor a(t) is governed by the Friedmann equation: 

^ a^ a-'^ 

where c is a constant. In what follows we will only be concerned with the case 7 = | corresponding 
to the so-called traceless perfect fluid (pure radiation). Furthermore, we assume A < 0. Equation 
([7|) admits a static (i.e. time independent solution) in which the values of the scale factor and 
the cosmological constant are related by: 

a{t) ~ oq ~ constant, A = Aq = (8) 

In the dynamical case, under the assumptions 7 = |, A < 0, the Friedmann equation ^ can be 
integrated explicitly — see e.g. [14]. Different types of solutions are obtained, depending on the 
relative value of A with respect to Ao as given in equation ([8]), where ao 7^ is now the value of 
the scale factor at some fiduciary time t = to 0. The relevant properties for the analysis of 
these solutions are summarised in the following proposition: 
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Proposition 1. For a FLRW cosmology with k — 1, ^ — ^ and A < 0, A 7^ Aq, the scale 
factor, a{t), is a smooth, non-vanishing and monotonically increasing function for t G [io,oo), 
with t = to > and oq — alto) > 0. Furthermore, 

ds 

Jto 

and one has the limits 

a — > 00, a/a^\J~-^\, ii/a— >— ^A. 
as i —>■ 00 . The pressure for these models is given by 

-4/4 

where po = p{t^). In particular, one has that p ^ as t ^ 00. 

The proof of this proposition fohows from direct inspection of the exphcit solutions — see e.g. 
[g, page 78. 

Remark 3. A similar type of result can be obtained for FLRW models with 7 = |, A < and 
k = -1, 0. Again, see 



4 The conformal field equations with matter 

The stability of the solutions to the Einstein equations described by the metric 5^1, corresponding 
to the line element ([5]) will be analysed in terms of a conformally related (unphysical) metric 
g^^. This strategy leads to consider the conformal Einstein field equations. The idea of vacuum 
conformal Einstein field equations expressed in terms of the Levi-Civita connection V of the 
metric and associated objects was originally introduced in [3J 1^1 IHI ■ The generalisation of 
these conformal equations to physical spacetimes containing matter was discussed in [lOj. More 
recently, a more general type of vacuum conformal equations — the extended conformal Einstein 
field equations — expressed in terms of a Weyl connection V has been introduced — see [11] . 



4.1 Conformal rescalings 

All throughout we assume that the two metrics g^,y and g^i, are conformally related to each other 
via 

ff^!^ = 6^3^!^, (9) 

where is a non-negative scalar field — the conformal factor. The Christoffel symbols T^p^ and 
r^/jy of the associated Levi-Civita connections V and V are related by 

f - F/, = S^^P^Tx, (10) 
where = and Sf^^^^ is the conformally invariant tensor 



4.2 Curvature tensors 

In a 4-dimensional spacetime the Schouten tensor, P^u^ of the connection V is defined by 

The Schouten tensor of the connection V is defined by a similar expression involving the physical 
Ricci tensor and scalar. The tensors P^^ and are related by 

P^. - P^. = V^T„ - T^T, + ig^.T.T". (11) 
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We can thus decompose the Riemann curvature tensor, R^vXp, of the connection V into its 
irreducible parts as 

= C'.xp + 2 [g'^ixPp]. - 9.[xPpf) , (12) 

where C^vXp denotes the conformally invariant Weyl tensor. 

As V is a Levi-Civita connection it satisfies the first and second Bianchi identities: 

R''wxp]=0, (13a) 
^[aR^'l.lXp] = 0. (13b) 

In our discussion of the conformal field equations with matter we will make use of the physical 
and unphysical Cotton- York tensors Y\pi, and Y\pi, given, respectively, by 

Yxp. = VxPp. - ypPxu, Yxp. = ^xPpu - ^pPxu. 
The tensor Iapi/ appears in the once contracted Bianchi identity 

V^C^Ap = n^.. (14) 
Finally, it is noticed that the twice contracted Bianchi identity takes the form 

"^"Ppu = VpP, (15) 

where P = g^^Pxt.- 

4.3 Frame and spinor formulations 

In what follows, consider a frame field {e^}, j = 0, . . . , 3 which is orthogonal with respect to the 
metric g^i,. By construction one has that 

gn^ei^ej" = r]ij, r]ij = diag(l, -1, -1, -1). (16) 

In order to discuss the extended conformal Einstein field equations, it will be convenient to 
regard, for the moment, the connection V only as a metric connection with respect to g^j^ — i.e. 
VAff/ji/ = 0. Under this assumption, the connection V could have torsion, and thus, it would not 
be a Levi-Civita connection. The connection coefficients, Fj'^j, of V with respect to the frame Ck 
are defined by the relation 

As a consequence of having a metric connection, the connection coeflScients satisfy 

Ti'^jVki + Ti'^irikj = 0. 
The torsion, Sj'^j, of the connection V is defined by 

If Sj'^j = so that the connection V is the unique Levi-Civita connection of g^y, the connection 
coefficients acquire the additional symmetry 

i 3 ~ j i- 

Related to the g'-orthonormal frame we will consider a normalised spinor dyad {6a}, A = 
0, 1, such that 

k 

SAA' = ekCT AA' 

where a'^AA' are the constant van der Waerden symbols. 
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In the sequel, a space spinor formalism will be introduced — see e.g. |25| . To this end, we 
consider a timelike spinor taa' which in terms of the dyad {Sa} can be expressed as 

AA' ^ A- A' , ^ A- A' 

In particular, one has the normalisation condition taa'T^^ — 2. The space spinor formalism 
allows to turn primed indices in spinorial expressions into unprimed ones by suitable contrac- 
tions with ta"^ — see [101 [TSl \T9\ [20] for more details. We simply recall that the space spinor 
decomposition of a spinor uaa' is given by 

UAA' = \UTAA' - t'^A'Uqa, (17) 

where 

_ ppi _ pi 

U = Upp'T , UAB^T (bUa)P'- 

4.4 The conformal field equations with tracefree matter 

In our subsequent discussion it will be convenient to distinguish between the geometric curvature 
r^iij — i.e. the expression of the curvature related to the connection coefficients Ti^j. — and 
the algebraic curvature R^uj — i.e. the decomposition of the curvature in terms of irreducible 
components given by equation (1121) . One has that 

r^,, ^ (r/,) - e, {t^\) - ^,„^ (r /'^ - r^™,) + r>„r^.™, - r/,„r ™, + s.-.r^^. 



Following [13], in the sequel it will be convenient to introduce the variables 

d^Uj = Q-^C''uj, (18a) 

= V,e, (18b) 

SEE i(V'=dfc + eP\.)- (18c) 

Furthermore, we also consider the following zero quantities — cfr. [lOJ: 

S,V,= (r/^-r/^)ez-[e„e,], (19a) 

=r^,, -i?^,,, (19b) 

Ahj = V,Pji ~ VjP,i - dkd^uj - Q^T.ji , (19c) 

Auj=Wkd''Hj~QT,,i, (19d) 

4 = 4- VfeO, (19e) 

S^J EE W^dJ + GPy - ST].,, - iO^Ty, (19f) 

Cfe = Vfc.s + d'Pki - ^e^d'Tik, (19g) 

C = X-GQs + Mkd'', (19h) 

where 

Tki = 6 '^Tki, Tijk = 6 '^Yijk- 



The interpretation of the zero quantities (|19ap - (|19d| is as follows: the zero quantity given 
by (|19ap measures the torsion of the connection V; that of (jl9b|) relates the expression of the 
curvature of V with its decomposition in terms of irreducible components. Equations (|19c[) and 
(|19dp measure the deviation from the fulfillment of the once contracted Bianchi identity. Finally, 
equations (|19ep . (I19f|) and p9gP bring into play the definitions (jl8bp and (|18cp and give rise to 
differential conditions for the fields Q, di and s. 

The conformal Einstein field equations with matter are then given by 

E^^efc = 0, = 0, A,.y = 0, Ahj = 0, (20a) 

4 = 0, 5,j = 0, a =0, C = 0. (20b) 
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These equations yield differential conditions for the frame coefBcients e^, the spin coefficients 
Ti^ k, the components of the Schouten tensor Pij, the rescaled Weyl tensor d'^uj, the conformal 
factor 0, the 1-form di, and the scalar s, respectively. As discussed in e.g. |T3], equation (jl9h[) 
has the role of a constraint which holds by virtue of the other conformal field equations if it is 
satisfied on some initial hypersurface. It is noticed that as the torsion, is being introduced 

as a zero quantity, it can be consistently set to zero in the geometric curvature appearing in the 
definition for the zero quantity E,^Hj — equation (jlQbll . 

Equations (j20ap - (|20bp need to be complemented with the energy-momentum conservation 
equation ([4]). Its particular details will depend on the matter model under consideration. 

Remark 4. Using a direct generalisation of the arguments presented in [5, 4J one can show that 
a solution to the conformal Einstein field equations with matter (|20ap - (|20bp and ^ give rise to 
a solution to the physical Einstein- matter system ©-(HI) — see also Theorem 3.1 in [6|. 

Remark 5. As a result of the conformal rescaling the conformal equations (|20ap - (|20bp have 
a built-in conformal freedom which needs to be specified in order to deduce suitable evolution 
equations for the conformal fields. Further gauge freedom in equations p0a[) - (l20bp is concerned 
with the partial specification of the frame ej. and the choice of coordinates. These will be specified 
by the choice of suitable gauge source functions. 



5 Perfect fluids in the context of the conformal approach 

In this section we present a discussion of the relativistic equations describing a perfect fluid which 
is geared towards our particular applications. 



5.1 The energy- momentum tensor and its transformation rules 

Given the spacetime {A4,g^^)^ the energy-momentum tensor for a perfect fluid with 4- velocity 
M*, pressure p, and density p has the form 

Tf^u = ip + p)Uf^u„ - pgi_,„. (21) 

In order to perform a discussion of the perfect fluid in the conformally rescaled (unphysical) 
spacetime one introduces unphysical versions of the physical fields. More precisely, one defines 

T^^ EE Q^^f^i,, EE eu^, p = Q^^p, p = Q^^p. 

Following the approach used in the discussion of geometric fields, we will work directly with the 
frame components Tij = ei^Cj-T^i/ and Ui = ei'^u^^ with respect to a y-orthonormal frame e^. 
Thus 

Tij (p + p)uiUj - pi]ij. 

We observe that g{u,u) = 1 implies that g{u,u) — 1. Now, using Ui = rjijU^ , — if-'uj, our 
signature convention implies 

— uq, u"" = —Ua, a = 1, 2, 3. 



A computation using the standard transformation rules for the covariant derivatives of con- 
formally rescaled metrics yields 

Consequently, the (physical) equation for the conservation of energy-momentum 

y'f,, = 0, 

implies an analogous equation 

V^T,, = 0, (22) 
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for the (unphysical) conformally rescaled spacetime only if the energy-momentum tensor is 
tracefree — see [10] . Notice that T = ff^ Tij = if and only if T = if'^ Tij = . A quick computation 
shows that for a perfect fluid the tracefreeness of the energy-momentum tensor implies p — 3p = 
— in other words 7 = |- Hence 

This class of perfect fluids is usually referred to as pure radiation. 

In the present article, our analysis will he restricted to the case of tracefree perfect fluids. The 
unphysical energy-momentum tensor for this class of perfect fluids reduces to 

T,j = ^u.Uj ~ \pr],j. (23) 

As a consequence of the definition of the 4- velocity u' it follows that 

rjijU^u-' — Ukv!^ ~ uqu^ + Uau"" — 1, (24a) 

Vfcu" = -^Vfe«^ (24b) 
Wo 

ViVfeu" - -^^i^^u'^ - L^^ua^^u- _ !^s/iu''\/ku-. (24c) 

Uq Uq Ug 

These identities will be used to rewrite the component and its derivatives in terms of the 
spatial components u'^ and their derivatives. This procedure will be central for the construction 
of a symmetric hyperbolic system for the matter variables. It is also noticed that equation (|24ap 
implies 

w'^Vfc {u,u') = 0. 

This expression shows that if UiU^ = 1 at some point in a fluid flow line, then UiU^ = 1 in the 
whole flow line. 

5.2 The energy conservation equation and the equations of motion 

A direct computation shows that the conservation equation ([22]) implies 

Zj = I [uju'-V.^p + pUjV.,u' + pu'V.,Uj) ~ ^Vjp = 0. (25) 

This equation can be split into components parallel and orthogonal to u': 

u'Zi = u'Vip + jpViu' = 0, (26a) 
j/Z, = Ipu'W^u, + luju'W.p - iV,p = 0, (26b) 

where 

7*i — ^ij UiUj. 

These equations are the conformal versions of the equation of energy conservation and the 
equations of motion — see e.g. |3J. It is noticed that equations (j26ap and (|26b[) can be combined 
to give 

VjP = Apu'ViUj - ^pujViu\ (27) 

This equation will be used in the sequel to eliminate the gradient of the unphysical density from 
certain expressions. 

5.3 A symmetric hyperbohc system for the fluid fields 

The equations of conservation of energy and motion will be used to construct a symmetric hy- 
perbolic system of evolution equations for the unphysical density p and the spatial components 
of the unphysical velocity m". The procedure used here follows the presentation given in f3]. In 
the sequel, it should be understood that, consequently with equation (|24ap . 

Uq — u'^ — \/l — UqW". 
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Substituting identity (|24b| into (|26a|) gives 

^ («°VoP + u'^Vap) + ^ (vaW^ - ^Voz.-^) = 0, (28) 

where the extra factor l/4p has been included to ensure symmetric hyperbohcity. Similarly, from 
equation (j26b[) one deduces 

ip Ap 

In order to obtain suitable evolution equations for the spatial components of the 4-velocity, we 
consider the combination 



or equivalently 



^'^ = U^.u'^ + ^^^^V.u^) + ^ f ^?7°'=Vfep - = 0. (29) 

\ u"mo / 4p y 

A direct inspection shows that: 

Lemma 1. Equations (j28p and (|29p constitute a symmetric hyperbolic system for the fields p 
and u"" . 

One also has that: 

Lemma 2. A solution {p, u"") to the evolution equations (j28[) and (j29[) implies a solution {p, uq, u°') 
to equation wif/i uq = \/l ~ UaU°-. 



Proof. We need to show Z-' = 0. The definition of uq implies that (j24bp and (|24c|) hold. Now, 
given a solution to ((28|) and ([29|) . the right hand side of (f28|) can be rewritten so as to yield 
Zj — 0. Substitution into the left hand side of ([29]) gives 

- = 0. 

Contracting with Ma and using (|24ap . as well as uq > 1 gives first Z'~' = and then = 0. Hence 
a solution to ((281) and ([29]) satisfies ((25)) . 

Remark 6. Let u-^^ denote the spinorial counterpart of the 4-velocity vector m''. The spinor 
u^^ can be split using the spinor as done in (fT7)) . This implies 

where cr'Xg denote the spatial Lifeld-van der Waerden symbols. It follows that (|25p implies a 
symmetric hyperbolic system for the spinorial components u and uab ■ The explicit form of these 
equations will not be required in our subsequent analysis. 

5.4 The Cotton- York tensor of a traceless perfect fluid spacetime 

The matter field quantities feedback into the geometric part of the conformal field equation 
through the physical Cotton- York tensor Yf^iuX- In what follows, the latter is expressed in terms of 
tensors, however the frame and spinor component versions are easily derived from these equations. 

For a tracefree energy momentum tensor the physical Schouten tensor is given by Pij — \Tij 
so that 
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Rewriting this expression in terms of unphysical quantities one obtains for = Q '^Yijk 
that 

Tijk = VfjTjjfe + TjjTjjfe + gfc[jTj];T'. (30) 

The last two terms in this expression are polynomial in p and the components u*. The first term, 
however, contains derivatives of and p that would enter the principal part of the Cotton- York 
and Bianchi equations. The fluid equations cannot be used to eliminate these derivatives. 

In order to get around this difficulty, we introduce new variables pk and uij and corresponding 
zero quantities qk and yij via 

qk = Pk-'^kP, yij = Uij ~ViUj. (31) 

Observe that if 5^ = and yij = 0, one then has that UijU^ = and yijU^ — 0, so that one can 
write 

WjO = Uia, yiQ = yia- 

Uq Uq 

Furthermore, from Uiju^ = 0, it also follows that 

kUi = VkUi H (2/jo-Uio)H [yia-Uia), (32a) 

Uq Uq Uq 

u^ViUkj = Ukjyi' ~ UkjUi^ ■ (32b) 

Finally, if = 0, yij = 0, then the first term of Tijk can be written as 

V[,;Tj]fe = |(V[jpWj]Mfe + pW[,Uj]Uk + pu[jW,]Uk) - lW[,pT]j]k, 

= |(P[iUj]Wfc + pU^ij^Uk + pUljUi]k) - ^PliVj]k- 

5.4.1 A symmetric hyperbolic system for pk and Uk"" 

The evolution equations for p and m° are derived from equation ([25]) . Taking derivatives of (|25|) 
and commuting them gives: 

= VfcZj = I {ujU'V^VkP + pUjV.VkU' + pu'V^SlkUj) - \VjVkP 

+^{ujU'Y,j}^Vip + pujrf./iu^ - pu'r^^^jUi + VuUju'Vip + VkpUj^iu' 
+Vkpu'yiUjUjVku'VtP + pVkUjVtu'' + pVku'ViUj) - iS^'^Vip, 



3 k j 

where all terms with at most one derivative of p or u'' have been gathered in Vkj ■ In view of this 
discussion, in the sequel we will consider the field equation for pk and Uij given by the following 
zero quantity: 

Zkj = I [uju'-ViPk + pUjViUk + pu'ViUkj) - ^VjPfe + VFfcj = 0, (34) 

with 

Wkj = ^{ujU^Y^kiPi + P'^j'Tkii'^'' - Pu^r^ijUi + Ukju'pi + pkUjU^ + pku'uij 
+UjUkPt + pukjUi^ + pUkUij) ~ \Y.^^jPi. 

From the equation Zij — one derives, in analogy to (|26b[) and (|26ap . that 

u^Zk, = u'^.pk + |pV,Mfe' + Xfe = 0, (35a) 
= ^pu'WiUkj + ^UjU^WiPk - I'^jPk + Xk] = 0, (35b) 
where (|32b| has been used and 

Xk = ^pu\ukjyi^ - UkjUi^) + u^Wkj, 
= -^puju\ukiyi - UkiUi) +j/Wki. 
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Finally, we rewrite (j35ap in the form 

KVoPfc + u'^V.pfe) + ^ (^V.wfc'^ - ^^ouk""^ +Xk=0 (36) 

where 

^ TTT^fe + 1~ (yiO - Mjo) H (Via - Uia) 



Similarly, the combination 



leads to the evolution equation 



u 



7 ^kl - 7 ^kl 

Uo 



u'^rUk'- + ^^-4^y^uA + ^ i^if^Wk - V^'^iPk ]+Xk'' = 0. (37) 



where is a combination of 

-7 -''^fei ~ 7 -^fe/ 



Wo 

and terms from expression p2ap . 

In analogy to Lemma [1] one can readily verify that: 

Lemma 3. If tjij — 0, then equations (j36p and (|37p constitute a symmetric hyperbolic system for 
the fields pk and Ufe". 

A similar argument to the one leading to Lemma [5] yields: 

Lemma 4. Let i/ij = 0. ^ solution {pkTUk"") to the evolution equations p6p and p7p implies a 
solution to equation 



5.4.2 The subsidiary equations for the fluid variables 

In this section we derive evolution equations for the zero quantities qi and yij . These subsidiary 
equations will be of relevance in the discussion of the propagation of the constraints — see Section 



Subtracting equation (1331) from equation (|34p gives 

Qkj = I {uju'V.iqk + pujV^yk' + pu'V^ykj) - ^"^jQk + Wkj - Vuj = 0. (38) 

Now, using substitutions like 

UijPk - ViUjVkP = Uijqk + VrjPk - Vijqk 

one can deduce that all individual terms in sums of W^j — V^j in equation psp contain at least 
one zero quantity. Repeating the discussion for the evolution equations for [pk^ Uij) with (qk, yij) 
one finds that 

{u^^oqk + u^'Vaqk) + ^ f Va2;fc" - — Voy^"^) + Ffe = 0, (39) 
lop^ Ap \ uq J 

u'V^Vk'^ + ^^^W.yk') + ^ f ^77°V,g, - rf'Wiqk) + Yk"" = 0, (40) 
u'^uq J Ap \u'^ J 

where all terms in Yj, and 1^." contain zero quantities. The evolution equations p9p and PP]) 
constitute a symmetric hyperbolic system for the independent components of qk and y^j . 
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5.4.3 Final remarks 

As a consequence of the analysis in the previous subsections one has that the components Tijk of 
the tensor T^i,\ with respect to the frame are polynomial expressions of the unknowns p, pj , 
Ui and Uij . If desired, the dependence with respect to pk can be eliminated using equation (|27|) . 
As long as p 7^ 0, the fields p, pj, Ui and Uij satisfy symmetric hyperbolic equation^j. Finally, it 
is noticed that if the fields p, pj, Ui and Uij are regular at the points where = 0, then Tijk is 
also regular — and consequently, also equations (|19cp and (|19d|) are formally regular. 

6 A symmetric hyperbolic reduction of the conformal field 
equations 

In the previous section it has been shown how the equations of motion for the fluid variables and 
their derivatives lead to a system of symmetric hyperbolic equations independently of geometric 
gauge considerations. The purpose of this section is to briefly discuss a reduction procedure for 
the geometric unknowns. Our treatment is inspired on the one given in [10], but it also combines 
ideas from |9l [19] . 

6.1 Gauge freedom 

As mentioned previously, the conformal Einstein field equations (|20ap and (|20bl) are endowed 
with three classes of gauge freedom: conformal, coordinate and frame. In what follows, we briefly 
discuss a procedure for fixing this freedom. 

6.1.1 Conformal gauge freedom 

As already mentioned, the conformal Einstein field equations (|20ap and (|20b|) admit certain 
freedom in the specification of the representative, g^^ , of the conformal class [g^u] which will be 
used as the unphysical metric — see e.g. |T3j and references therein for more details. 

Assuming for a moment that one has a solution to the conformal Einstein field equations with 
matter, it follows then by contraction of indices in equation ([TT]) together with the tracefreeness 
of the energy momentum tensor that 

As discussed in, e.g. [T5], this equation can always be solved locally so that the condition 

P = Pfc^- = -l, (41) 
holds. This condition fixes the conformal freedom in the equations (|20ap and (j20bp . 

6.1.2 Coordinate and gauge freedom 

In order to fix the coordinate and frame gauge freedom, we make use of the notion of gauge source 
functions — see [TJ I10[ I13j . The gauge source freedoms will be chosen so as to render symmetric 
hyperbolic evolution equations for the geometric unknowns. 

In what follows we encounter equations of the form 

2V[,M,]... = iV[,,]... 

where the dots denote an arbitrary set of indices — cfr. equations (|19ap - (|19cp . The spinorial 
equivalents of the above equation are given by 

V^(^-MV)... (42) 

^An alternative evolution system which is valid even if p = can be obtained using the ideas of |22 |. 
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and its complex conjugate. Now, an equation of the form 



(43) 



is well known to imply a symmetric hyperbolic evolution system for the independent components 
of Maa' — see e.g. [7]. Note, however, that equation (j42l) contains no information about the skew 
term 

Va[A'M^B']- - ieA'B'V^^'Mcc'... = F..., (44) 

which can be specified arbitrarily. Thus, by adding (|44|) with a convenient choice of a gauge 
source function, F..., to (|42| one obtains an equation of the form of (|43|). from where a symmetric 
hyperbolic system can be extracted — see e.g. [71 113|. 

The previous discussion will be implemented in the field equations (|19ap . (|19bp and (|19cp . 
These equations provide differential conditions for the fields ef , Pij and Ti^ k- Let e^^, , Paa'BB' 
and Taa'^^ cc denote the spinorial counterparts of these fields. As a consequence of the metric- 
ity of the connection, instead of working with Taa'^^ cc we will consider a spinorial field 
^AA' BC such that 

BB' -r B ^ B' , B' ^ B 

i AA' cc — AA' C^C +^A'A C'^C ■ 

For convenience, define the gauge source functions 

F^^W^^'e^,, (45a) 

F(^BC) = y^-^'rAA'BC, (45b) 

Fbb' = V^^'Paa-bb- = Vbb'P, (45c) 

where the second equality in the definition of Fbb' follows from the twice contracted Bianchi 
identity for the unphysical connection V. Motivated by their value in the reference solution (the 
conformal FLRW solution) the gauge source functions will be fixed by the conditions 

^" = 0, F(BC)=0, Fbb-=Q- (46) 

Notice, in particular, that the last condition is consistent with the conformal gauge condition 
(|4T|l . As discussed in — see also [10] — a particular choice of the coordinate and frame gauge 
functions F^ and F(^bc) fixes the coordinates and fram^l- 

The spacetimes to be considered in the present analysis have the topology of K x S'^. Given 
an initial manifold S for the spacetime, then there is a diffeomorphism $ : 5 — ?> which allows 
to pull-back coordinates from to S. These coordinates on the initial manifold S will be used 
as the initial value of the spatial part of the spacetime coordinates. The time coordinate will 
be set initially to zero. The initial value of the frame is set by choosing on S some arbitrary 
orthonormal spatial frame Ca (with respect to the 3-metric of S) . The cq vector is set to coincide 
initially with the (spacetime) normal to S. 

6.2 The evolution equations 

The hyperbolic reduction of the matter variables has already been discussed in Sections 15.31 
and I5.4.T1 In what concerns the evolution equations for the geometric variables, we follow the 
procedure indicated in [TUj. This consists of a rewriting the spinorial version of the conformal 
field equations (|20ap - (|20bl) in terms of space spinors so that the resulting equations contain only 
unprimed indices. In order to encompass the full information of the field equations, one has to 
include into the set of equations their Hermitian conjugates. If the fields and equations are then 
decomposed into their irreducible parts, then the equations split in a natural way into symmetric 
hyperbolic evolution and constraint equations. This procedure is straightforward, but involves 
lengthy computations, most of which can now be implemented in a computer algebra system like 
the suite xAct for Mathematical. 

^The gauge source functions _F'' and F^bC) implyi respectively, wave equations for the coordinates and a 
semi-linear equation for the frame components. These equations can be solved locally. 
•^See www.xAct.es. 
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The required evolution equations have aheady been deduced in [lOj. Their detailed form will 
not be required here. Instead we present a summary of their key structural properties. In what 
follows let 

V = (9, dAB,s, e^QjTABCD, Pabcd) , 

<t> = 4>ABCD, 
Q = {P,U^AB)) , 

■0 ^ {pab,uab(cd)) , 

where only the independent irreducible components of the spinors are taken into account. In 
terms of these objects, the evolution equations have the form 

90^^ + Aj^^{v)cry = B^^}(v)v + M[„](^>, 0, £>, -0), (47a) 

{y2E + Af^j {v)) doct> + Af^j {v)crj> = {v)4> + M[^] {v, xjj, g), (47b) 

Af^] (i^, Q)daQ + (r, Q)crg = B[^] {v)g, (47c) 

A°^] (w, Q)doip + Af^j {v, Q)crip = iv)ip + M[^] {v, cj), q). (47d) 

In equations (|47al) - (j47d[) . E denotes the 5x5 unit matrix, while A-^j, A-p^j, A||^j denote smooth 

symmetric matrix-valued functions of their respective arguments. In particular, A|^j (0) = and 

A|^j , A|^j are positive definite if p > 0. In addition, , , -B[g] and JB[^] are smooth matrix- 
valued functions of v. Finally, Af[^,], iVfj^j and A/j^j are non-linear vector- valued functions of 
their respective arguments. These functions are smooth if u'^ ^ 0. 

For convenience of the discussion, in the sequel, we define 

w = (Re(u),Im(u),Re((/)),Im((/)),Re(£.),Im(£.),Re(t/>),Im(i/j)). 

Remark 7. The deduction of the evolution equations (|47ap - (|47dp assumes the choice of gauge 
source functions given in (|46l) . 

6.3 Propagation of the constraints 

An analysis of the so-called subsidiary equations describing the propagation of the zero quantities 

Si'j, S'^Hj, ^lijj ^lij: ^fel Ck, C> (48) 

in terms of which the geometric part of the conformal field equations — see equations (j20ap - 
(I20b[) — is expressed has been given in [TU]. This lengthy analysis is succinctly summarised in the 
following lemma: 

Lemma 5. If the unphysical energy-momentum tensor, Tij, satisfies 

V'Ty - 

and the expressions 

V'T,,fc, (ie^T-d^.fc, + VSTkij + eV.Tfc,,) - d^e^\^^Tun. (49) 

can he rewritten in terms of matter zero quantities, then the geometric zero quantities in (|48p 
satisfy a subsidiary system which is symmetric hyperbolic and homogeneous in the zero quantities. 

A lengthy computation assuming the form for the energy-momentum tensor given by equation 
and taking into account the expression ([50)) for the rescaled Cotton- York tensor, shows that 
the expressions (HIJ)) in Lemma [S] can be rewritten as a homogeneous expressions of the matter 
zero quantities q^, yij, Zkj and Qkj defined by equations (|3ip . ([34]) and ([38)), respectively. The 
analogue of Lemma [5) for the matter zero quantities is given by Lemmas [2] and |4] 

The purpose of the analysis of the propagation of the constraints is to establish the follow- 
ing reduction theorem which follows directly from the symmetric hyperbolicity of the subsidiary 
systems and their homogeneity with respect to the zero quantities. 
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Theorem 1. A smooth solution w of the propagation equations (|47ap - (j47b[) which satisfies the 
constraint equations on a spacelike hypersurface S defines in the domain of dependence of S a 
solution to the conformal Einstein field equations with matter model given by a traceless perfect 
fluid. 

7 The traceless perfect fluid FLRW cosmology as a solution 
to the conformal Einstein fleld equations 

The purpose of the present section is to cast the traceless perfect fluid FLRW cosmology with 
A < in a form in which its character as a solution to the conformal Einstein field equations with 
matter becomes manifest. 

7.1 The FLRW cosmology on the Einstein cyhnder 

One of the most important properties characterising FLRW cosmologies is their conformal flatness. 
This shows that as in the case of the Minkowski, de Sitter and anti-de Sitter spacetimes, these 
solutions admit a conformal representation in which the unphysical spacetime {M,g^u) is given 
by the so-called Einstein cylinder (or Einstein cosmos). 

The Einstein cosmos is given by the manifold A^^ = R x S'^ with a metric given by the line 
element 

gg = dr^ - dcr^ (50) 

where, again, da^ is the standard line element of S'^. The manifold S'^ will be coordinatised in the 
way indicated in Subsection [221 A f;-orthonormal frame Ck can be defined on Mg by completing 
the frame {ci, C2, C3} discussed also in Subsection 12.21 with the vector eo = dr. Setting, for 
convenience, cq = 9r, one can write — Ck-Cs with efe- = 5k- — the components of Ck with 
respect to the basis Cg. 

In order to relate the FLRW line element ([6]) with that of the Einstein cosmos, equation ([50|. 
one introduces the change of coordinate 

ds 

a{s)' 

This naturally leads to the following choice of conformal factor: 

e(r) ^ l/a(r) ^ l/a(t(r)), (51) 

so that gg = Q^g^, where g^ is given by equation ([5]). For the class of FLRW cosmologies 
covered by Proposition [T] one has that G) — > as t — > 00. Moreover, there exists a (finite) positive 
constant Tqo such that &{too) — 0. Notice also that r = for < = to. 

A direct computation using the line element (|50l) . the frame and the conformal factor ([5T|) 
gives the following expressions for the unknowns of the conformal field equations: 

4- = 4-, fP^^^eoilkV^^, P^j^5-^5j° -\ri,j, d,jki = 0, (52a) 
e = a-\ dk - -a-^a'Sk", s = ^a'^a'^ - la'^a" - la-\ (52b) 
p = poQq, III = Si°, pk = 0, Uij = 0. (52c) 

where ' denotes differentiation with respect to r, and oq and po are the values of the scale factor 
and the physical pressure at the initial time r = 0. Notice that the unphysical density for 
this model is constant. Spinorial versions of the above expressions can be readily obtained by 
contraction with the constant spacetime Infeld-van der Waerden symbols a^AA', or their space 
spinor version a°'AB- The explicit expressions will not be required in our subsequent analysis. 
Following the notation of Lemma [7l we collect the independent spinorial components of the fields 
in (|52ap - (|52cp in a vectorial unknown which we denote by w. 
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A direct computation using the expressions (|52a|) - (|52c|) shows that, for this solution, the gauge 
source functions F^, F(^ab) and Faa' as defined by (|45ap - (|45cp are given by 



F'' = 0, F^AB) = 0, Faa' = 0. 
This computation justifies the choice of gauge source functions made in 

Recalhng that d/dr — ad/dt, and using the hmits given in Proposition [U it foUows that 

dk — y^— |A, J ^ as T Too- (53) 

Accordingly, the expressions given in (|52al) . (|52b|) and (|52cp define a smooth solution to the 
conformal Einstein field equations p0a|) - (j20bp for r G [0,Too]- In fact, this solution extends, at 
least locally, beyond t = Too- This can be easily seen to be the case by using the values of 
the solution ()52a|) - (|52c|) as the initial value for a Cauchy problem on the slice t — Too- The 
symmetric hyperbolicity of the evolution equations implies that the solution to this initial value 
problem exists for r G [too,Too + S) for some 6 > 0. From the expression for dk in (j52b[) and 
Proposition [1] it follows that do < at Too- Thus, by continuity, S can be chosen such that 9 < 
on (too, Too + S) - In summary we have: 

Lemma 6. There exists S > such that the expressions (|52a[) - (|52cp give rise to a solution to the 
evolution equations implied by the conformal Einstein field equations (|20a[l - (|20b[l on [0,Too + S). 
Furthermore O < m (too,Too + <^). 

Initial data for a FLRW cosmology 

The expressions in ()52ap - (|52c|) naturally induce an initial data set for the conformal Einstein field 
equations which we denote by Wq. Notice that there is no need for performing a pull-back in this 
construction as the fields in (|52al) - (|52cl) are all scalars. 

7.2 Structure of the conformal boundary 

The structure of the conformal boundary for the solution to the conformal Einstein field equations 
described by (I52ap - (|52cp follows directly by inspection. 

By construction at t = Too one has that 8 = 0. From the limits ([53|l — see also equation 
(jl9hp — one has that: 

^d'^ = VfcOV'^e = -iA > 0, at T = Too, 
so that the future conformal boundary ^+ = {p E M.g | = 0} is spacelike. 

8 Existence and stability results 

The purpose of the present section is to provide our main results. These concern the global 
existence of solutions to conformal Einstein field equations with matter source given by a perfect 
fluid in the case A < 0, 7 = | which can be regarded as non-linear perturbations of the reference 
FLRW solution described by Proposition [TJ We also provide results concerning the structure of 
the conformal boundary for these solutions. Altogether these results show the non-linear stability 
towards the future of the reference FLRW cosmological model. 

8.1 An Ansatz for the solution 

We will consider solutions to the evolution equations (|47ap - (j47d|) of the form w = w + w , where 
w as defined in Section [7. 11 and w describes a non- linear perturbation from the reference solution 
w. The fields in w are interpreted as the pull-back of the original fields on Mg under a cylinder 
map. In what follows let Wq — Wq + Wq be an initial data set for the system of evolution 
equations (|47ap - (j47d|) prescribed on the initial manifold S — as discussed in Subsection 18.21 It 



17 



will be assumed that Wq satisfies the conformal constraint equations. The vector Wq is to be 
interpreted as the pull-back of the smooth map relating and the initial manifold S. 

A direct inspection gives rise to the following result: 

Lemma 7. For w sufficiently close to and as long as p > and vP ^ 0, the equations (j47ap - 
(j47dp imply a symmetric hyperbolic evolution system 

3 

A^(w + ih) ■ drW + ^ A-{w + w) ■ Cr{w) + B{t, X, w, CgW, w) ■ w = 0, (54) 
1=1 

for the independent components ofw. The matrix valued functions A-, B are smooth Junctions 
of their arguments. Furthermore, the entries of the matrix A^{w) are hounded from below by 
l/%/2. Finally, w = Q is a solution of equation (|54p . 

8.2 Constructing initial data for the conformal evolution equations 

In the sequel, it will be assumed that one has a solution {S, hap, K^p, p,il'^) to the (physical) 
A < Einstein-perfect fluid constraint equations 

P + - Kapk"'^ ^2{X^ p), (55a) 

b'^k^p-bpk^y, (55b) 

with S having the topology of and the perfect fluid satisfying a barotropic equation of state 
with 7 = |- In equations (|55ap - (j55bp bp and f denote the Levi-Civita covariant derivative and 
the Ricci scalar of the intrinsic 3-metric hap of S. Kap is a symmetric 3-diniensional tensor 
corresponding to the extrinsic curvature of 5 with respect to the 5- unit normal n^. Furthermore, 
jl = h'^n'^Tfj^^, while jp corresponds to the pull-back to S of j\ = n^^h\^T^^. For a perfect fluid 
with a tracefree energy-momentum tensor, a direct computation gives that 

jl = ■|p(4m|| - 1), in = Ipup^, 

where -uy = u'^h^ and up corresponds to the pull-back of hfi'^Uu to S. In particular, if on S one 
has that and are aligned — as in the case of the FLRW cosmologies — then jl = p and 
jp = 0. In general, however, we will consider perfect fluid configurations for which and u'^ are 
not aligned. 

Using a generalisation of the procedure for vacuum spacetimes described in, say, [6] one can 
construct a solution to the conformal constraint equations implied on S by equations (|20ap - (j20bl) . 
Following the notation introduced in Section [6.21 we denote the independent components of such 
a solution by Wq. 

8.3 The main result 

In what follows, given m e N, let || • ||m denote the Sobolev-like norm on the space C°°(§^R^) of 
smooth valued functions on for some non-negative integer N — see e.g. [8j [18] for precise 
definitions. Furthermore, let iJ"'(§^,R^) be the Hilbert space obtained as the completion of the 
space C•°°(§^R^) in the norm || • \ Using the cylinder map between Ms and M, one can 
apply the norm || • ||„i to evaluate the norm of functions on the unphysical initial hypersurface 
S. Furthermore, the vector w ~ w{t, x) can be regarded as a function of r which takes values in 

Our main result is the following theorem: 

Theorem 2. Suppose m > 4. LetS denote a 3-dimensional manifold diffeomorphic to S'^, and let 
Wq = Wq + Wo be initial data for the conformal evolution equations (j47ap - (|47dp constructed from 
some physical initial data set, (S,hap, Kap, p,u°'), for the Einstein field equations with A < 
and matter source given by a traceless perfect fluid (7 = 4/3). There exists e > such that if 
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Il'i'ollm < £ then the initial data setwo determines a unique solution, w, to the evolution equations 
(|47ap - (|47dp which exists on [0,r*] with > Too- The solution w is of class C™^^([0,t*] x S'^) is 
such that: 

(i) it determines, in turn, a C*™^^ solution to the A < conformal Einstein field equations, 
equations p0ap - (j20b|) and p5p - (|34p . with matter given by a traceless fluid on [0,t,] x S^; 

(ii) there exists a function Tj^ — Tj^{x), a; G S"^, such that < T+(a;) < and 

e > 0, on 7W = {(T,a;) e M X §3 I < r < r+(a?)}, 

6 = 0, dkd^ = -\X<Q on = {(r+(a;), a;) e R x | a; e S^}. 

(Hi) one obtains a (7™^^ solution (M,gfj_^, p,u^), to the A < Einstein-perfect fluid field equa- 
tions with 7 = 4/3 which is future geodesically complete for which ^+ as defined above 
represents conformal future infinity; 

(iv) given a sequence of initial data iUq"^ such that \\w^^\\m < £ and \\wl^'^\\m — as n — ^ 
c», then for the corresponding solutions ti;^"^ (with minimum existence time T*j one has 
I I'U'''"'' I |m — > uniformly in t £ [0, r*]. 

Remark 8. The above theorem, and in particular part (iv), amounts to a non-linear stability 
result for the A < 0, 7 = | FLRW cosmological models, in the sense that sufBciently small 
perturbations of data for the FLRW solution give rise to (future) global solutions to the Einstein 
field equations with the same asymptotic structure as the reference solution. 

Remark 9. Note that no consideration of the vorticity of the radiation fluid was required for 
the derivation. The vorticity of the fluid can be calculated from the components of Uij. 

Proof. Existence, uniqueness and the smoothness of the solutions to equations (|47ap - (|47dp follow 
from the properties of the equation (j54p provided in Lemma [7] and an extension of the general 
existence and stability Theorem by Kato ^6J provided in [8J — see also [18]. In particular, if 
e is sufficiently small, one obtains a common existence time > Too for all initial data with 
ll'ii'ollm < £■ Part (iv) follows from the same result. 

Now, the Reduction Theorem, Theorem [T] ensures that if the conformal constraint equations 
are satisfied on S, one obtains a C""~^ solution on [0,t*] x S"^ to the A < conformal Einstein 
perfect fluid equations with 7 = | — this shows part (i). 

In order to show part (ii), one observes that because < in (too,Too + S), then if e is 
sufficiently small one has that 8 < at, say, t = Too + S/2. As 9 > at t = 0, then there is a 
for which B = 0. By reducing, if necessary e one has that such t is unique, and hence, the function 
T+(a;) is weU defined. If follows from and (|T91il) that 9 = implies dkd*" = VfeOV'^e = 

— iA > 0. Hence 8 = 0, respectively t = T+(a;), defines a regular spacelike hypersurface 

For part (iii) one notices that a solution of the vacuum conformal Einstein field equations 
implies a solution to the vacuum Einstein field equations — see e.g. [6nl0j. If w''' denotes the dual 
cobasis of the frame e^, (w'^, e^) = S'^i, then the unphysical metric is given hy g — (E) lu-' ■ If 
e is sufficiently small, one has that det(g) 7^ on [0,t*] as det(g) 7^ 0. Thus g = Q^^g is well 
defined on Ai. The matter fields p and are defined via the formulae in (PT|) . An adaptation 
of this argument to our setting gives the desired result. Geodesic completeness follows from an 
analysis of the geodesic equations and standard perturbative arguments for ordinary equations 
given that the background solution is future geodesically complete. 

□ 
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